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ABSTRACT 

We consider generically the problem of a Gamma-ray burst (GRB) fireball with an additional energy 
injection, either in the form of a Poynting-flux-dominated outflow or a kinetic-energy-dominated matter 
shell injected after the burst. Generally, a total injection energy comparable to that of the impulsive 
energy in the initial fireball is required to make a detectable signature in the afterglow lightcurves. 
When this criterion is met in the case of Poynting-flux-dominated injection, this leads to a gradual 
achromatic bump appearing in the otherwise power-law afterglow lightcurve. Alternatively, in the case 
when the injection is kinetic-energy-dominated, the results depend on whether the collision between the 
rear (injected) and the forward shell is mild or violent. If the relative velocity between the colliding shells 
does not exceed a critical value defined by their energy ratio, the collision is mild, and the injection may be 
analogous to the Poynting-flux injection case. Otherwise, the injection is violent, and an additional pair of 
shocks will form at the discontinuity between two colliding shells, so that there are altogether three shock- 
heated regions from which the emission contributes to the final lightcurves. We describe the shell-merging 
process in detail including collision and relaxation by taking into account the dynamical evolution and 
the emission from the various shocks involved. Assuming synchrotron emission, we calculate afterglow 
lightcurves in the X-ray, optical and radio bands for the various cases. The injection signatures due 
to violent matter-dominated collisions are abrupt and complicated, due to the emission from any of 
the three emitting regions, and depending on the injection parameters and the observed energy bands. 
This differs from the gradual bump signature found in the Poynting-flux injection case. In both the 
Poynting-flux-dominated and the kinetic-energy-dominated cases, the energetics of the fireball as well 
as the absolute afterglow flux level after the injection are boosted with respect to the one without post- 
burst injection. Identifying the different injection signatures from future early afterglow observations 
may provide diagnostics about the nature of the fireball and of the central engine. 

Subject headings: gamma rays: bursts - radiation mechanisms: non-thermal - shock waves - stars: 
magnetic fields 



1. INTRODUCTION 

Classical Gamma-ray burst (GRB) afterglow models 
generally invoke an "impulsively" injected (possibly col- 
limated) fireball running into an ambient interstellar 
medium (ISM) or into a pre-burst stellar wind. In princi- 
ple, after the initial fireball starts to decelerate, additional 
injection of energy is still possible. There are at least three 
motivations to study the post-burst injection possibility 
more carefully. 1: Additional injection is a natural expec- 
tation of the internal shock GRB model, which predicts 
that some slow shells will trail the outermost fast shell 
and run into it when the latter is decelerated (e.g. Rees & 
Meszaros 1998; Kumar & Piran 2000). In some cases, the 
central engine may also inject fast shells at a later time. 
Later injection is also naturally expected in certain types 
of the central engine models (e.g. Dai & Lu 1998; Zhang 
& Meszaros 2001a). 2: A GRB fireball can in principle 
be dominated by the kinetic energy of the baryons or by 
a Poynting-flux component. Any clue about the fireball 
nature is presently lacking. As discussed below, an injec- 
tion signature in the GRB lightcurve may directly provide 
diagnostics about the nature of the injection as well as in- 
formation about the central engine. 3: "Bump" features 
have been observed in several GRB afterglows (e.g. GRB 
970228, GRB 970508, GRB 980326, GRB 000103C, etc.). 



and various interpretations have been proposed, e.g., re- 
freshed shocks (Panaitescu, Meszaros & Rees 1998), su- 
pernova components (Bloom et al. 1999; Reichart 1999; 
Galama et al. 2000), dust cchos (Esin & Blandford 2000), 
and gravitational micro-lensing events (Garnavich, Loeb 
& Stanek 2000). It would be of great interest to study the 
injection feature in more detail to find distinct properties 
to be differentiated from other interpretations. Such injec- 
tions are likelier to occur in the early afterglow phase. The 
planned future broad-band GRB mission Swift will have 
the ability of recording broadband early afterglow signals 
from many GRBs, and thus bring unique opportunities to 
study the injection features. 

Depending on the different types of the central en- 
gine, the post-injection could consist of, e.g., some kinetic- 
energy (i.e. baryon) dominated shells, or a Poynting-flux- 
dominated wind (Usov 1994; Meszaros & Rees 1997b). 
Here we present a generic study on various injection cases. 
Generally speaking, an injection energy comparable to the 
initial impulsive energy is required to make the injection 
noticeable. We therefore study a system containing an im- 
pulsive shell which is already heated during the shell-ISM 
interaction, and which is collecting material from the ISM, 
and in the meantime also receives a large enough injection 
energy from a continuous Poynting wind or a trailing mat- 
ter shell. We first discuss in §2 the simplest case, in which 
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the fireball continuously receives pure energy with negligi- 
ble baryon loading, very likely in the form of Poynting flux. 
Since no reverse shock is expected in such a case, the injec- 
tion signature is solely from the forward shock emission, 
through the change of the blastwave global dynamics. In 
§3, we discuss the more complicated matter shell injection 
case. In real situations, shell collisions may occur more 
than once. Here we only investigate the detailed physics of 
one such post-collision. The same analysis, when applied 
several times, is generic in delineating the more compli- 
cated multi-collision cases. We first analyze the general 
hydrodynamics of the three-shell (impulsive shell, injec- 
tive shell, and ISM) interaction (§3.1). After reviewing 
the general shell evolution history, we solve numerically 
for the condition of a "violent injection" , in which an ad- 
ditional pair of strong shocks form at the discontinuity of 
the two colliding shells (§3.2). We then outline in §3.3 
the hydrodynamics of the three-shell-interaction process 
in detail by dividing the whole process into five stages 
(§3.3.1-§3.3.5), namely (1) pre-coUision; (2) mild injection; 
(3) violent injection; (4) relaxation; (5) post-relaxation. In 
§3.4, we discuss synchrotron radiation from all the possible 
shocked regions in various stages, and present the injection 
signatures in the broadband lightcurves. Our findings are 
summarized and discussed in §4. 

2. INJECTION FROM A POYNTING-FLUX-DOMINATED 

FLOW 

The conventional definition of a Poynting-flux- 
dominated flow is L±^^^js[ / Lp <C 1, where £±,-y.Ar is the 
luminosity of the electron-position pair, radiation and 
nuclcon components, and Lp is the luminosity of the 
Poynting-flux (i.e. magnetic field) dominated outflow 
component (e.g. Usov 1999). In the context of GRBs, 
discussions of a Foynting-dominated flow as compared to 
a kinetic-energy-dominated flow are mainly focused on 
the paradigm of an approximately baryon-free Poynting 
flow, which are not as conducive to coUisionless shock 
formation, so in this case GRBs may be powered by the 
non-linear breakdown of large-amplitude electromagnetic 
waves (Usov 1994; Blackman & Yi 1998; Lyutikov & 
Blackman 2001; for a review, see Usov 1999). However, 
in both scenarios, coUisionless external forward shocks are 
believed to form in the afterglow phase (Usov 1999). Al- 
though the baryon free nature at the GRB prompt phase 
is uncertain (e.g. Spruit, Daigne & Drenkhahn 2001), in 
the later injection phase the bayron loading could be in 
principle much lower and the injection could be in the form 
of an almost bayron-free Poynting-flux-dominatcd wind. 
In such a case, no coUisionless reverse shock is expected to 
propagate into the wind (Kennel & Coroniti 1984, Usov 
1999; essentially, because the sound speed in the ejecta 
is close to the speed of light), and the injection energy is 
used to increase the total energetics of the fireball. 

Considering an adiabatic relativistic hot shell which is 
collecting material as it expands into the ISM and which 
receives an energy input from the central engine via a 
Poynting-flux-dominated flow, the differential energy con- 
servation equation in a fixed frame is (i[r(Moc^ -I- M2C^ -f 
U)] — dM2C^ +dEinj, where Mq, and M2 are the masses of 
the impulsive shell and the swept-up ISM, respectively^, 

'^The subscript 2 is used to match the notation in the kinetic-energ; 



£'inj is the received injection energy, and F is the bulk 
Lorentz factor of the blastwave. The injection wind can be 
approximated as being cold, so it does not provide any in- 
ternal energy, and the total internal energy in the fireball is 
provided by the shock heated ISM. Noticing equation (||), 
one can usually approximate [/ ~ (F — 1)M2C^ as a work- 
ing assumption to perform a first-order treatment of the 
problem. Strictly speaking, such an approximation is only 
valid when the post-shock fluids are considered uniform, 
i.e., with constant density, temperature, and Lorentz fac- 
tor. More detailed self-consistent modelings require tak- 
ing into account the Blandford-McKee (1976) self-similar 
profiles and tracking the post-shock internal energy more 
closely. With the assumption, the integrated energy con- 
servation equation then reads 

Eini + (To - T)Moc' = (F^ - l)M2c'. (1) 

This is analogous to the combination of equations (1) and 
(2) in Zhang & Meszaros (2001a). When E'inj is negli- 
gible and noticing that iSimp — FqMoc^, equation (|^) is 
the standard equation for the impulsive adiabatic blast- 
wave evolution (e.g. Chiang & Dermer 1999; Huang, Dai 
& Lu 1999), which can be reduced to the familiar form 
Ei^p = (47r/3)i?3„TOpc2F2 (so that F cx R'^/^ cx T-^/», 
where R is the blastwave radius in the fixed frame, and 
T is the observer time) when Fq ^ F ^ 1. In the pres- 
ence of ii'inj, we can see that the condition for changing 
the blastwave dynamics is that Ei^j exceed (Fq — F)Afoc^, 
or essentially i?inj 

A Poynting-flux-dominated wind is usually continuous, 
with an intrinsic luminosity law which may be defined 
through a temporal index q within a certain time regime, 
e.g. L{T) (X T^, where T is the intrinsic time of the cen- 
tral engine (which is the same as the observer's time after 
the cosmological time dilation correction). A detailed de- 
scription of such an injection scenario, especially in the 
context of magnetar central engine models, was discussed 
in a previous paper (Zhang & Meszaros 2001a). Here we 
summarize some of the conclusions. When the blastwave 
dynamics is dominated by the injection (i.e. iJinj^-Eimp, 
which only occurs when q > —1), the bulk Lorentz factor 
evolves according to F cx i?-(2-«)/2(2+9) ^ y-(2-g)/8_ 
most problems, usually q — during the injection phase, 
so that F oc c>c T~^/^. There are three relevant 

time scales in the problem, i.e., the blastwave deceleration 
time, To; the critical time when the continuous- inject ion 
energy component dominates the impulsive-injection en- 
ergy component, Tc] and the characteristic time durung 
which the central engine produces a substantial energy 
imput (i.e. q > —1), T. The condition for the injec- 
tion signature to show up in the afterglow lightcurves is 
T > max(rc,To). For < To < T, the dynamics is 
defined by the continuous-injection as soon as the after- 
glow is set up, and we define it as an immediate injec- 
tion (Fig. la). The lightcurves will be flat from To to T, 
and will steepen after T. For Tq < Tc < T, the sig- 
nature will not show up until Tc, and we define it as a 
delayed injection (Fig. lb). There are two temporal in- 
dex changes: the lightcurves will fiatten around Tc, and 
resume the original steepness around T. The injection- 
signature is achromatic due to the change of the blast- 

'-dominated case discussed in §3, e.g., Figure 2 and equation (|L2|). 
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Fk;. 1. Injection lightcurvcs from a Poynting-flux-dominatcd flow ejected by a highly magnetized mihisecond pulsar 
(magnetar) central engine. Solid and dashed lines are the cases with or without a continuous energy injection. Shown 
are the bulk Lorentz factor of the blastwave, the total energy in the fireball, as well as the lightcurves in various bands 
(X-ray: v = 10^^ Hz; optical: f = 10^^ Hz; radio: u = 10^° Hz). Parameters adopted: .Bimp = 10^^ ergs, n = 1 cm~^, 
Tq = 300, Eg — 0.1, e_B ~ 0.01, and p = 2.5. Different pulsar parameters lead to the injection being either "immediate" or 
"delayed", (a) Immediate injection: Bp = 10^^ G and Pq = 1.5 ms; (b) delayed injection: Bp = lO^^G and Pq = 1.5 ms. 



wave dynamics. Since no reverse shock is expected, only 
emission from the forward shock is responsible for the af- 
terglow emission. Assuming synchrotron emission of the 
relativistic electrons, the temporal index at the injection 
phase is a = (1 — q/2)(3 + q + \ (which is /3 -f 1 for 5 = 0), 
as compared to the standard case of a = (3/2)/3, where /3 
is the spectral index in a certain band of interest. As an 
example, we again discuss the Poynting-flux injection case 
in the magnetar model. In this particular case, the injec- 
tion luminosity law is L{T) = Lpm.oi^ + T/Tem)~'^, where 
Lem,o = ^ 1-0 X lO^^erg S-1B2 i^PoT-s-Re and 

= 3c^I/B^R^nl = 2.05x103 s h^Bpl^Pl.^R-^'' , and 
Bp = lO^^cmSp^is, Pq = 10^''sPo,-3 are the pulsar dipolar 
field strength at the pole, and the initial period at birth, 
I = lO^^g cm^/45 and R = lO^cmiJg are the rotation in- 
ertia and the radius of the pulsar, respectively. Figure 1 
shows two example injection lightcurves for both the im- 
mediate and the delayed injection cases in three energy 
bands. A general feature is that the injection signature 
is rather smooth, mainly due to the smooth varying lumi- 
nosity law L{T). Dai & Lu (2001) and Wang & Dai (2001) 
obtained a similar result recently. 

3. INJECTION FROM A KINETIC-ENERGY-DOMINATED 

SHELL 

The previous Poynting-dominated injection is, for the 
present purpose, the simplest case, in the sense that the in- 
jection does not result in new emission sites. The signature 
in the lightcurves is still produced by the shocked ISM, al- 
beit through a changed global dynamics of the blastwave. 
However, if the injection energy is dominated by the ki- 
netic energy of the baryons, collisionless shocks may in 
principle form at the discontinuity of any two adjacent 



shells under certain conditions. The emission sites during 
the injection process therefore coidd be multiple, and the 
injection signatures, due to the joint contributions from 
all these sites, should be much more complicated. In this 
scx:tion, we investigate in some detail the shell collision 
process, as well as the emission signatures. 

3.1. Hydrodynamics 

We consider a collision in the afterglow phase, in which 
case the impulsive shell is already heated up by the passage 
of the reverse shock. Assuming that the electron equipar- 
tition factor Cg is not close to unity, the impulsive shell 
will not cool at the time of injection, even if the elec- 
trons cool rapidly. A generic configuration of the problem 
is summarized in Figure 2, which is predicated upon the 
fulfillment of certain conditions (§3.2). There are six re- 
gions of interest, namely: (1) unshocked ISM; (2) forward- 
shocked ISM; (3) imshocked hot impulsive (leading) shell; 
(4) forward-shocked impulsive shell; (5) reverse-shocked 
injective (trailing) shell; (6) unshocked injective shell. We 
denote through Mi, rii, and pi the mass, baryon number 
density, internal energy density, and the pressure of the re- 
gion "i" in its own rest frame, respectively; 7^ and are 
the bulk Lorentz factor and the dimensionless velocity of 
region "i" measured in the fixed frame, respectively; and 

7ij - hi/lj +7jhi)/^, for » 1 (2) 

and (3ij the bulk Lorentz factor and the dimensionless ve- 
locity of region "i" measiired in the rest frame of region 
"j", respectively. In the hot regions (2, 3, 4 and 5), we as- 
sume a relativistic equation-of-state, i.e., Pi = {j — l)ei = 
ei/3, where 7 = 4/3. Two contact discontinuities separate 
the three parts, i.e., the ISM, the impulsive shell (which we 
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Fig. 2. — A schematic (not to scale) description of the violent injection configuration from a kinetic-energy-dominated 
shell. The thick solid line indicates the pressure, or the internal energy density, in each region. The thin solid line indicates 
the co-moving mass density in each region. 



will consistently call shell 3 in the following discussions), 
and the injective shell (which we call shell 6 in the follow- 
ing discussions). The pressure equilibrium and velocity 
continuity at the contact discontinuities yield 



where 
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There are three shocks that accelerate particles and con- 
tribute to the emission, i.e., (A) a first forward shock prop- 
agating into the ISM, (B) a second forward shock propa- 
gating into the impulsive shell, and (C) a reverse shock 
propagating into the injective shell. The continuity of the 
energy, momentum, and particle flux densities in the rest 
frame of each shock leads to the following jump conditions 
(Blandford & McKee 1976; Sari & Piran 1995; Kumar & 
Piran 2000) 



— = (72 - I)r7ipc2, !^ = 472 + 3, (5) 
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2 (1 + 364/63) (3 -I- 64/63) 
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(64/63)(l + 364/63) 
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(7) 
(8) 



= (756 - l)TOpC , = 4756 + 3 . 

715 ne 

Notice that the jump conditions at the second forward 
shock (eqs. |^ and j^]) are different from the others, 
mainly because the region in front of the shock (region 
3) is hot. With (|) and (|), a combination of (|) and @ 
yields 



64 "6 f (755) £^0,6 £"(756) . /, R 

— = ; — — ~ • • mm 1, 

63 ni F(73) Eo,3 4 V R^^q 



(9) 



£(7) = (7- 1) (47 + 3) ~ 472, 7>1. (10) 



The second half of the equation (^) will be explained fur- 
ther below. Given 73 and Tg, the unknown 74 can be solved 
for by replacing (^) into (g) with the use of (||). 

3.2. Shell evolution & condition for violent collision 

Equations (H)-(^) imply that strong shocks form be- 
tween the contact discontinuity of the injective shell and 
the impulsive shell (e.g. Blandford & McKee 1976), as in- 
dicated in Figure 2. A physical solution for the non- linear 
equations (^-(§) requires 64 > 63, or 73 < 74 < 75, which 
demands a minimum relative Lorentz factor between the 
two colliding shells given the parameters of both shells 
(cq.[^). This is qualitatively in agreement with Kumar & 
Piran (2000), who demand a minimum enthalpy density 
ratio for a given 735 (in their Fig. 3), which is equivalent to 
our case in which we demand a minimum 735 for a given 
Eq,q/Eo^3. For a shock to form, the injective shell should 
move super-sonically with respect to the impulsive shell, 
i.e., by noticing the sound speed — e/\/3 for the shell 
3 (in the relativistic hmit 7 = 4/3), 736 > ^3/2 ~ 1.22. 
Usually this condition is automatically satisfied given the 
strong shock forming condition 64 > 63. In this subsection, 
we examine the criterion for 64 > 63, which is the condi- 
tion for the configuration in Figure 2 being realized. We 
first summarize the evolution of a matter fireball shell (e.g. 
Meszaros, Laguna & Rees 1993; Piran, Shemi & Narayan 
1993; Kobayashi, Piran & Sari 1999; Piran 1999), and 
then show how a shell collision process is well-defined by 
the initial conditions of the shells. 

The injection of a matter shell from the central en- 
gine can be quantified by three independent parameters, 
i.e., the energy of the shell Eq, the mass of the shell 



5 



Table 1 



Shell evolution: Evolution of the thickness (A in the fixed frame, A' in the comoving frame) and the comoving density n 
of a shell as a function of the fireball radius R. no = Mo/[mp(47r/3)Ao] is the initial baryon number density, rii is the 
density of the unshocked ISM, and M2 is the mass of the shocked ISM. 



condition 


radius regime 


A' 


A 


n 


n/ni 




Ao <R<Rc 


R 


Ao 


Mo 
^7r_R^ rrip 




M0/M2 


Rs < Rd 


Rc ^ R ^ Rs 


Rc 


Ao 


( ^" ) ( ^) 


= no(^)^ro 


{Mo/M2)iR/Rc) = iEo/M2c')iR/Rs) 


(thin shell) 


Rs < R < Rd 


R/To 


R/Tl 


( Mo A p 


(Mo/M2)To = S0/M2C2 




R> Rd 


R/T 


R/T^ 


( M„ \ p 

\ ^TrR^nip J 


= no(^)^r 


(Mo/M2)r = (£;o/M2c2)(ro/r) 


Rs > Rd 


Rc < -R < Rd 


Rc 


Ao 


( Mo \ ( R\ 
y^TTR^mpJ \RoJ 


= no(^)^r 


{Mo/M.2){R/Rc) = {EolM2^){RIRs) 


(thick shell) 


R > Rd 


R/T 


i?/r2 


( Mo A p 


(Mo/M2)T = (Eo/M2C^)iTo/T) 



Mq (thus the dimensionless entropy, or the maximum 
Lorentz factor of the shell is rj = Tq ~ Eq/Mqc^), and 
the initial thickness of the shell Ao which is essentially 
determined by the timescale of the central engine activ- 
ity, i.e., Aq — cST. The following parameters can be 
consequently defined: the initial baryon number density 
of the shell, which is no = Mo/[mp(47r/3)Ao]; the ra- 
dius at which the fireball stops accelerating and enters 
the coasting phase, which is Rc — FqAo; and the ra- 
dius at which the fireball starts to spread with its local 
sound speed due to the velocity difference within the fire- 
ball matters, which is Rs — FqAq. There is another pa- 
rameter, i.e., the deceleration radius Rd, which usually 
appears in the problem. However, this radius not only 
depends on the intrinsic parameters of the shell (e.g. to- 
tal energy), but also on the properties (e.g. density and 
velocity) of the medium that stops the shell. For the inter- 
action between the impulsive shell (shell 3) and the ISM, 
the deceleration radius is roughly the radius where the 
mass of the ISM collected by the fireball is 1/70,3 of the 
shell mass, i.e., Rd.3 ~ R-y.s = (3i?o,3/47rnimpC^77|)^/'^ = 
1.2 X lO^^cm £^d,3^5inj;^/^(7o,3/300)-2/3^ where ni is the 

ISM baryon number density in unit of Icm"'^, -©0,3,51 = 
-Bo,3/10^^ergs (hereafter Xn = X/10", except those sub- 
scripts that denote different regions, such as ni). For the 
interaction between the injective shell (shell 6) and the 
impulsive shell (shell 3) , the deceleration radius is around 
the radius where the injective shell catches up with the 
impulsive shell. Let us denote AT36 as the time inter- 
val between ejecting the two adjacent shells (as compared 
to ST, which is the timescale of the injection process), 
the "catching-up" radius is then Rd^ ~ 27|c(Ar36)/(l — 
7I/7I). For 76 ^ 73, Rdfi is mainly defined by the 
Lorentz factor of the slow, leading shell, i.e., i?d,6 — 

27|c(AT36) = 2.1 X IQi^cm E^oll^^{AT3e/20sy/^n;^^\ 
The condition for the collision to occur in the afterglow 
phase is then Rd,e > Rd.3, or approximately Ar3g > 

2.3s El^l^^n~^^^{jo,3/300)-^/^. The deceleration radius 
of a shell can in principle be cither smaller or larger than 
the spreading radius, Rs- For Rs > Rd (a "thick shell", 
e.g. Sari & Piran 1995), the shell has not started to spread 
when the deceleration (or the collision) occurs. Alter- 
natively, for a "thin shell" (with Rs < Rd), the shell is 
spreading with its sound speed in the comoving frame due 



to the velocity difference within the shell as the collision 
occurs. In this sense, the initial "thin shell" is actually 
thick during the collision. Table 1 summarizes the evolu- 
tion of the shell widths (both in the fixed frame. A, and in 
the comoving frame. A'), shell baryon number density n, 
and the ratio n/ni in various regimes. In the last column 
(i.e., n/ni), for the rear shell 6, we have assumed that the 
radius of the shell 6 front is comparable to that of the shell 
3 front, i.e., we assume that A3 ^ R, which is generally 
the case as long as 73 » 1. We note that expressing n/ni 
in terms of the mass (or energy) ratio of the shell and the 
collected ISM (M2) allows a clearer physical understand- 
ing for the problem. At the coUsion point Rd.e, it is conve- 
nient to show that ng/ni = {M0.6/M2) ■ min(i?/_Rc.6j 7o,6) 
(thick shells for the former, and thin shells for the lat- 
ter). Noticing that 70^/0,0 — £'inj — Eq,q, and that for 
7o,3 ^ 73 3> 1 (i.e., the coUision occurs long after the 
afterglow's setting-up but well before the non-relativistic 
phase, which is usually the case), one has ©(73) — 47I, 
and 7IM2 ~ i?imp = Eo^3 (see eq.Q and the discussions 
thereafter). This justifies the second part of the equation 

The condition for 64 > 63 can be investigated numeri- 
cally. Figure 3 shows the minimum 735 required for 64 > 63 
for different energy ratios of the two shells, i?o,6/£'0,3- No- 
tice that the wiggles on the curves are due to the crudeness 
of the numerical solutions. A general trend is that a larger 
relative velocity between the two shells is required for a 
less energetic trailing shell. The critical line depends on 
the status of the shell 6, i.e., whether it is expanding (thin 
shell) or not (thick shell) when the collision occurs. This 
is least demanding for an expanding (thin) shell 6, which 
is plotted as the solid curve. For a thick shell 6, the larger 
the ratio Rs/R, the more demanding the condition is. This 
is because for the same energy ratio £'o,6/-^,3 7 a thicker 
shell tends to get a smaller 64/63 ratio (eq.|^). The dot- 
ted and the dashed curves in Figure 3 are for a shell with 
Rs.6 5 and 10 times the collision radius, Rd,6, respectively. 
Another comment is that, as long as in the relativistic 
phase lasts, only the energy ratio and the relative velocity 
are relevant to the problem. The absolute values of the 
shell energies and Lorentz factors do not enter the solu- 
tion. Notice that the supersonic condition defines another 
critical line for the shock forming condition (Fig. 3). The 
case for 64 > 63 not being satisfied is less clear, but in 
this case no strong shock is expected. We regard such an 
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injection case as "mild" (with weak or non-shock), and we 
expect that the injective materials are effectively attached 
onto the forward shell to increase the total energy budget, 
which is analogous to the Poynting-flux injection case. 

In summary, the input parameters in the injection prob- 
lem include: (i) Ei^p = Eq^s = 7o,3^o,3C^; (ii) 7o,3 
(or A/o,3, which is only important around the on-set of 
the afterglow); (iii) ni; (iv) Einj = Eq^ = jofiMofiC^; 
(v) 7o,6 (or Mofi); (vi) Ao,6 = cSTq, which determines 
Rs,6 = 7o,6^o,6; (vii) ATae (or equivalently Rdfi). The 
comparison between i?s,6 and Rd,6 then decides whether 
the trailing shell is "thin" or "thick" when the collision oc- 
curs. Given a set of injection parameters (i) - (vii), at any 
radius R, one can first check the shock forming condition 
in Figure 3. If the condition is not satisfied (below the crit- 
ical line), the injection is a "mild" one. If the condition is 
satisfied (above the critical line in Fig. 3), the shell-merging 
process is violent. 

3.3. Shell merging process 

Though an injection of matter shell could be mild or vi- 
olent throughout, in some cases the injection could be mild 
in the beginning, and turn violent at larger radii when 73 
drops so that 735 meets the critical value. In this sub- 
section, we present an example which includes both cases 
during the injection. There are four characteristic radii in 
the problem: (1) the collision radius, Rcoi = Rd,&] (2) the 
radius where the persistent shocks form (Fig. 3), Rsh] (3) 
the radius where the injection process ceases, i?A; (For a 
mild injection, this is the radius where all the material in 
the shell 6 catches up with and "attaches" onto shell 3; for 
a violent injection, this is where the reverse shock crosses 
the rear shell and the merging process enters the relax- 
ation stage.) (4) the radius where 64/63 drops to unity 
and the relaxation finishes, Rj. The whole shell-merging 
process may be then be (at most) divided into five stages 
(Fig. 4a): (1) pre-coUision, R < Rcoi', (2) mild injection, 
Rcoi < R < Rsh] (3) violent injection, Rsh < R < Ra- 
(4) relaxation, i?A < R < Rf, and (5) post-relaxation, 
R> Rf. Below we discuss each stage in some detail. 

3.3.1. Pre- collision stage 

Let us consider a central engine that ejects an impulsive 
firebalP (shell 3) with Ei^p = i?o,3 = 7o,3-Mo,3C^- Some 
time (AT36) later, it ejects another (effective) shell (shell 
6) with Einj = Eofi = 70,6-^^0, 6C^- The forward shell starts 
to decelerate at i?d,3 due to the interaction with the ISM, 
and at a certain time before collision, it has collected a 
mass M2 worth of ISM, and has a Lorentz factor of 73. 
The rear shell moves faster, and we assume that before 
collision it is essentially undecelerated, and is still propa- 
gating with 7o,6- The whole system has two independent 
sub-systems. The deceleration of the forward shell is quan- 
tified by (cf. eq.0 but without the injection term) 

(70,3 - 73)Mo,3 = (73' - 1)^2 , (11) 

which is the standard adiabatic afterglow evolution. The 
rear shell is not noticeable at this stage. 

^The impulsive fireball could be in principle highly erratic, as indie 
collide and merge into one big shell, and we regard this as equivalent 



3.3.2. Mild injection stage 

Starting from Rcoi — Rdfi — '^liic^T^a), the shell 
6 catches up with the shell 3. At Rcoi < R < Rsh, 
the injection is not energetic enough to excite shocks. 
The injection is analogous to the case of a Poynting-flux- 
dominated flow, i.e., mainly it increases the blastwave 
energetics. Assuming that at a certain radius (time) a 
mass of M5 from the total injection shell Mg^o reaches 
the shell 3 and its energy is added to the total blastwave 
energy (here the subscript 5 is only for easy cornparison 
with the following notations, eqs. ||, in and ll, but 
does not have the physical meaning illustrated in Fig. 2), 
one can write down the global energy conservation of the 
system. The total energy of the two shells and the (col- 
lected) ISM if there were no interaction among them is 
(7o,6Afo,6 + 7o,3-^o,3 + Ai"2)c^; that with the interactions is 
joAMo^e-M5y+-f3{M5c^+Mo^3C^ + U2+3+5)- These two 
should be equal. Noticing that U2+3+5 — (73 — 1)M2C^, 
one has 

(70.6 - 73)M5 + Mo,3(7o,3 - 73) = hi - l)M2 . (12) 

This is very analogous to (||). For 76,0 » 73, the first term 
on l.h.s is approximately 76,0-^^5 1 which is just the energy 
injected into the system. 

During this mild injection phase, the mass in the shell 
3 gradually increases while that in the shell 6 (un-coUided 
part) gradually decreases. This ends when the 64 > 63 
condition is satisfied at Rsh and a pair of shocks (B and 
C) propagate into the two colliding shells, or, in a less 
energetic situation, when the whole shell 6 is "attached" 
onto the shell 3 and the injection process ceases. For 
the former case, we define Esh.3 = £^0,3 -I- M^c"^ and 
Esh,e — Eofi — Mc^c^ , both of which will be used later 
in the stage 4 (§3.3.4). 

3.3.3. Violent injection stage 

If the shell 6 is energetic enough or fast enough (see 
Fig. 3), one then has the configuration shown in Figure 2. 
The two new shocks (B and C) propagate into the two 
colliding shells, respectively. In principle, one can define 
this stage to end when either of the two shocks crosses the 
corresponding shell. However, the reverse shock always 
crosses the rear shell first, and we define the radius where 
this happens as i?A- The violent injection stage is then 
defined as Rsh < R < Ra- 

To describe the physical process at this stage, we need 
to quantify the dynamical evolution of the three shocks, 
or equivalently, the temporal evolution of 73 = 72 and 
75 = 74. Strictly speaking, 73 and 74 should be solved 
self-consistently, and one needs one more independent re- 
lation between 73 and 74 besides (^ and (joj, in order 
to get a self-consistent solution. However, this is usu- 
ally not necessary at this stage. Again, we write down 
the global energy conservation of the whole system, i.e., 

70,6Mo,6c2 + 70,3^0,36^ + M2C^ = 70,6(^0,6 - M^^^ + 
74(M5c2-HM4c2 + C/4-H5) + 73(M2c2-HMo,3c2-M4c2-hi72+3). 

Noticing C/2+3 - (73-1)^26^, U^+b - {lz&-l)M^c'^ , this 
gives a time-dependent equation 

(70.6 - 74756)^4 ~ (74 - 73)A^4 

+ (70,3 - 73)Mo,3 = (73' - ^)M2 ■ (13) 

!d by spiky GRB lightcurves. In the afterglow phase, many mini-shells 
the central engine efi^ectively ejecting one single shell. 



0.01 0.1 1 10 100 

Fig. 3. — The condition for a violent injection. The three curves are the critical lines for 64 > 63 being satisfied, for 
different cases of the rear shell (thin or thick at the collision). The super-sonic condition (the horizontal line) is also 
plotted. The space region above both the critical line and the horizontal line is where a violent collision as indicated in 
Fig. 2 occurs. 



Comparing to equations (|T^) and (0), we find that the 
first term in those equations (the already injected energy) 
is now replaced by the first two terms in ( [l3|) . The global 
dynamics (73 evolution) is modified only when the sum 
of the first two terms is comparable to the third term 
which is essentially i?imp/c^ — 70,3-^^0,3- This is impos- 
sible, since the energy already injected into the system 
(i.e. 70,6-^5) is now partially consumed in maintaining the 
kinetic energy of the injected material [i.e., —74756-^50^, 
since for 7o,6^i73 the term (70^6 ~ 74756) is nearly zero], 
and in raising the Lorentz factor of the newly shocked re- 
gions [i.e. the term —(74 — 73)1/4]. Unless i?inj ^ Ei^p, 
which is rather unlikely in GRBs, the sum of the first two 
terms is negligibly small, and ([l3| ) is essentially reduced 
to (^1]). The simple scaling law of 73 for the impulsive 
injection case is then applicable at this stage, e.g., for an 
isotropic adiabatic fireball running into a uniform ISM, 
one has 73 cx cx T~*/^. 

The unknown quantity 74, of interest here, can then 
be solved for by using (|^) and (^). The temporal evolu- 
tion of both 734 and 755, as well as the random Lorentz 
factor in the region 4, 7^, are thus well quantified. Al- 
though generally this is done numerically, a rough an- 
alytic estimate is possible. When the solution of 74 is 
physical, one can approximate 75 ^ 74 ^ 73, as long as 
76/73 > 1- In this hmit, 756 76/274, 743 ~ 74/273, 
and one has 74 cx (n6/ni)^/*7g^^. For our assumption 
of 76 ^ 76,0 =const, we have 74 oc n^^^ oc R^^/'^ for a 
spreading shell (cx R~^/^ for a non-spreading shelP). We 
can further estimate 64/63 cx 7I3 cx Ug 
(X R^). Since 63 cx 7I cx R^^ , we have 64 cx rig^^ cx R^^/'^ 
(or cx R^^). With (0), we also have 714/713 cx (64/63)^/'^ cx 



/^3/4 j^Qj. ^ Noticing that the shell 3 is spreading and 
that 7i3 (X i?~^73 cx one then has 714 cx R~^^/^ 

(or cx i?~^/^). The "random" Lorentz factor in the re- 
gion 4 (at the shock B) which defines the radiation from 
the region (notice that it is different from 743) is then 

73 = e4/7747npc2 cx nJ^R^ cx R^'^ (or cx i?5/2). Finally, 
the random Lorentz factor at the shock C is essentially 

7c = 756, which scales as cx n^^^^ cx i?'^/'* (or cx R^^'^). 
These very crude estimates generally agree with the nu- 
merical results (see Fig. 4a). 

Finally we estimate i?A- A self-consistent calculation 
of the shock crossing time (for both shells) requires hy- 
drodynamically tracking the speeds of the shocked and 
un-shocked shells. Very roughly, the shock crossing time 
(e.g., for the rear shell 6)_may be estimated (Sari & Piran 
1995) iA, = [Ag/c(/36 - M ■ (1 -76n6/75n5) ~ (A6/c)7| 
(for 76 ^ 75 = 74 ^ 73), where x denotes the average 
value of the quantity x during the passing process. This 
estimation is only good for a not very long passing time, so 
that the variables do not change much during the cross- 
ing period. This is the case for the rear shell, but not 
for the forward shell. For example, for a thin rear shell, 
A ~ R/je, and tAg ~ (-R/c)(74/76)^- This is much less 
than the dynamical time of the shell propagation in the 
fixed frame, R/c, as long as 74 ^ 7g. For a thick rear 
shell, the crossing time is even shorter since the shell is 
not spreading. For the forward shell, however, one has 
^ (^^/c), which is very long. Therefore it is evident 
that tAe is always shorter than iAs , and we use it to de- 
fine the end of the collision stage. The radius Ra is then 
defined as Rcoi + ct As , and we denote the corresponding 
parameters at this radius with a subscript "A", e.g., 6a, 3, 
eA,4, 7A,3, 7A,4, etc. These quantities are recorded from 



^Hereafter the scaling for an spreading shell is not bracketed, while that for a non-spreading shell is. 
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the code and are inputs to the dynamical problem of the 
next stage (see below). 

3.3.4. Relaxation stage 

Beyond i?A, the reverse shock passes across the rear 
shell, and the process enters the dynamical relaxation 
stage. After this point, there is no further significant en- 
ergy input into the region (4+5). Since 64 > 63 is still 
satisfied at Ra, the forward shock B still exists and prop- 
agates into the forward shell as long as the region (4-1-5) 
is separated from the region (2-1-3) super-sonically. The 
region (4-1-5) then expands adiabatically, since there is 
practically no heat flow across the shock. The internal 
energy density in the region (4-1-5) decreases according to 

the adiabatic law, i.e., 64 (x U4+5/V' oc {V')^^ oc R~'^'y^^^ 
for 7 = 4/3. As a result, 64/63 as well as 743 drops rapidly 
towards unity. During this process, 73 increases and 74 
drops, until the whole region (2-|-3-f4-t-5) reaches the same 
velocity, and contacts sub-sonically. The shock (B) dissap- 
pears and the relaxation stage ends, at the radius defined 
by Rf. 

During the relaxation stage, the shock conditions (g|), 
(||) and (Q) still hold. Again one needs to quantify the dy- 
namical evolution of 73 and 74. Unlike the previous stage, 
73 no longer evolves according to the simple oc R^^/^ form. 
Let us again write down the global energy conservation of 
the whole system, i.e., 70, 6-^0, ec^ + 7o 3-^0, ac^ + M2C^ — 

74(Mo,6c2 + MiC^ + Ui+5) + 73(M3,oC^ - M4c2 + M2C2 + 

C/2+3). Noticing that the shock heating energy produced 
at (A) is only consumed in the region (2-f3), i.e., U2+3 — 
(73 - l)M2c2, this gives 



(7o,( 



74)^^0,6 
+ (70,3 - 



- (74 - 73)^^ - 
73)^^0,3 = (7^ 



74C^4+5 
- 1)M2 . 



(14) 



This is even more complicated than (|l3|). As the relax- 
ation goes on, 74 gets closer to 73 and drops rapidly, 
thus both the second and the third terms decline. The first 
term, which is essentially the injected energy, starts to in- 
fluence the blastwave dynamics, just like the Poynting-flux 
injection or the mild injection (stage 2) cases (cf. eqs.|]l| 
and H). 

The relaxation is a complicated process. Strictly speak- 
ing, one needs to solve (|lj) and (||) to find a self-consistent 
solution for 73 and 74. Here we perform a simpler but ad- 
equate treatment. Since 73 is boosted up relative to the 
non-injection case, we assume that 



7I cx 63 cx R 



-3+5 



(15) 



where 6 is an unknown positive index. At Rf, the in- 
ternal energy density of the whole fireball is then 63J = 
eA,3{RA/ Rf)^~^ ■ Since beyond Rf the fireball evolves as 
if the initial total energy is i?o,6 + -E'0,3 = E^hfi + Esh.,3, 
then 63J can be also expressed in an adiabatic expan- 
sion form, but as if the energy density at i?A is a fac- 
tor (1 -I- Esh,6/Esh,3) times larger, i.e., 63J = (1 + 
Eshs/ Esh,3)eAARA/ Rff ■ This gives 



In 1 



(16) 



The last task is to determine Rf. In the relax- 
ation stage, 63 can be expressed as = eA,3(^A/^)'^^'', 
according to the assumption made. On the other 
hand, 64 evolves according to the adiabatic law, i.e., 
64 = eA 4(-Ra/-R)^(74/7a 4)'*^'^. One then has 64/63 = 

(eA,4/eA!3)(i?A/i?)(^+')/^(7A.3/7A,4)^/3(74/73)'/'. At 
R = Rf, 64/63 = 1 and 74/73 = 1 by definition. Noticing 
again equation (E5h, the relaxation stage ends at 



Rf — i?A 



eA,4 
eA,3 



1/3 / X 4/9 

^7A£ 

7A,4 



sh.6 



sh,3 



-1/9 



Since the parameters at i?A are in principle known 
(§3.3.3), equations (|5l)-(|l3), together with (|) and (|), 
adequately describe the dynamics in the relaxation stage. 

Some crude estimates match the numerical results. 
With 64 oc R-^it^^, 63 oc i?-(3+<5) and 7I3 oc 64/63 
(eq.[||), one can get, roughly, 74 oc a very steep 

declination (Fig. 4a). As long as the shock (B) still exists, 
new electrons are accelerated, and strong emission from 
the region heated by the shock (B) continues throughout 
the relaxation stage. At the rear end of the shell (the re- 
gion 5), since there is no longer a reverse shock (C), no 
new electrons are accelerated. However, the already accel- 
erated electrons are still hot and radiate synchrotron emis- 
sion. The random Lorentz factor in the region 5, which we 
still denote as 7^ (although the shock (C) no longer exist), 

also drops rapidly. Since 65 = 64 oc (V')^^ oc i?~*74^^, 
and 715 oc {V')^^ oc _R^'^74, one has 7^, oc [V')^^^^ oc 

i?~^7y^ oc R~'^ (Fig. 4a). Since 64/63 no longer 3> 1, there 
is no simple scaling-law estimation of 7^ at this stage. 

3.3.5. Post-relaxation stage 

After the relaxation ends, one has 73 = 74, the origi- 
nal 6^4+5 has died out, and the new thermal energy pro- 
duced at the shock (A) is consumed in the whole shell 
(2+3+4+5), i.e., f/2+3+4+5 - (73 - l)M2c2. The global 
energy conservation now reads 7o,6-Wo,6C^ + 7o,3-Mo,3C^ + 
M2C^ = 73(Afo,6C^ + Ms^oc^ + C/2+3+4+5), which is 



(70,6 - 73)^^0,6 + (70,3 - 73)^/0,: 



(73^ 



1)M2 . (18) 



This can be also obtained by removing the second and 
third terms in ( [T^ ) , as required by the physical relaxation 
process. Equation (|l8|) is a standard adiabatic expansion 
equation with a total energy of jo.eMofiC^ + 70,3^0,36^ = 
Einj + £'imp- At this stage, the only shock that accelerates 
new electrons is the shock (A) , and the emission from it is 
well described by the standard afterglow theory. 

Although shocks (B) and (C) no longer exist, their 
already accelerated electrons are still present and emit- 
ting. The random Lorentz factors in both regions (4 
and 5), which we still denote as 7^ and 7^,, evolve as 

oc (F 



/\-7+l 



OC 



R-^jI'^ OC i?-3/2 (Fig.4a). 



3.4. Synchrotron radiation and injection lightcurves 

With the shell-merging dynamics quantified above, we 
can calculate the injection signatures in the afterglow 
lightcurves. We have developed a code to calculate the dy- 
namics of the shell merging process as well as the emission 
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Fig. 4. — An exapmle of the kinctic-energy-dominatcd shell injection case. Parameters: for the impulsive shell 3, total 
energy E'o.s = 10^^ ergs, dimensionless entropy 70,3 = 300; for the injective shell 6, £^0,6 = 5 x 10^" ergs, and 70,6 = 500. 
The central engine activity time for the rear shell is 6Tq = 100 s, the time interval between ejecting the two shells is 
AT36 = 10 s. ISM baryon number density ni = 1 cm~^, electron spectral index p = 2.5. Equipartition parameters 
€e = 0.1, and cb = 0.02. The source redshift z = 1. (a) Evolution of various Lorentz factors with blastwave radius R. 'jj^, 
7^, and 7^, are the random Lorentz factors of the emission sites heated by the shocks (A), (B), and (C), respectively. The 
curves are solid when the shocks are on, and are dotted when they are off. 73, 74, and 75 are the bulk Lorentz factors 
of the regions (2+3), (4+5), and 6, respectively (cf. Fig. 2). Four critical radii, i.e., Rcoi, Rsh, Ra and Rj, separate the 
process into 5 stages, as marked in the plot, (b) The injection lightcurve in the optical band {v = 10^^ Hz). The emission 
components from all the three sites as well as the total flux are plotted. 




T (second) 



10-' - 




T (second) 



Fig. 5. — Broadband injection signatures (X-ray, optical, and radio) for the injection of a kinetic-energy-dominated shell. 

Solid lines arc lightcurvcs with injection, compared to the dashed line case without injection, (a) A "violent" injection 
involving formation of new shocks. All the parameters are the same as adopted in Figure 4. (b) A mild injection. Injective 
shell parameters are £^0,6 = 1 x 10'^^ ergs, 70.6 = 200. The central engine activity time for the rear shell is STq = 20 s, 
the time interval between ejecting the two shells is AT36 = 20 s. Other parameters are the same as in Fig.4 and Fig. 5a. 
Notice that although the total energy for the shell 6 is higher in the present case, the relative Lorentz factor between the 
two shells is much smaller, and the injection is mild. 
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from various shocks. Since the normal afterglow emission 
is generally well described by synchrotron radiation, we 
have neglected the high energy spectral components, such 
as the electron inverse Compton and the emission from the 
protons. Although these components could be important 
in the MeV to GeV range (see Zhang & Meszaros 2001b 
for a detailed discussion), usually synchrotron emission is 
adequate to describe the low energy afterglow (X-ray, op- 
tical and radio), at least in the early afterglow phase when 
injection may be expected. Nonetheless, the electron cool- 
ing due to the synchrotron self-Compton is incorporated 
through the Y parameter (see, e.g. Panaitescu & Kumar 
2000; Sari & Esin 2001; Zhang & Meszaros 2001b). 

The random thermal energy in a certain region is an es- 
sential parameter in calculating the synchrotron spectrum 
and the lightcurves. In stages 1 and 2, the only emitting 
site is the region heated by the shock (A). Starting from 
the stage 3, the emission sites increase to a total number 
of three, corresponding to the three shocks (although the 
shocks (B) and (C) die out in the stages 5 and 4, respec- 
tively). When the shocks are on, the random Lorentz fac- 
tors at the shocks (A) and (C) are essentially the Lorentz 
factor of the shocks, i.e., 7^ =73, and 7^ = 756- For the 
shock (B), since the unshocked region (region 3) is also hot, 
the random Lorentz factor is not simply 743 , but is rather 
7b = e^/nimpC? . When the shocks (B) and (C) die out, 
the random Lorentz factors in the relevant regions are still 
defined as 73 = e4/n4mpC^ and 7^, = e^/n^mpC?, which 
evolve adiabatically as discussed in §3.3.4 and §3.3.5. To 
calculate the synchrotron emission, we assume uniform 
values of the equipartition factors, i.e., cb and Cg, as well 
as the same electron power law index p in the three shock- 
heated regions. This is a crude approximation, but it is 
the simplest case which allows a straightforward compar- 
ison of the relative importance of the three shocks. The 
magnetic field in the region i is then Bj = (STresCi)^/^. As- 
suming that the injection factor of the electrons is unity, 
the characteristic injected electron Lorentz factor is then 
7m,x = ee[(p- 2)/ ip- l)](mp/me)7x, where 7^ is the ran- 
dom Lorentz factor at the region X (A, B or C), mp/rrig 
is the proton-electron mass ratio. Cooling frequencies are 
determined from the cooling electron energy, which is de- 
fined by the equality between the comoving time and the 
electron radiative cooling time. For the regions where the 
shocks are on and keep accelerating electrons, the syn- 
chrotron emission spectrum is calculated according to the 
four-segment broken power law description (Sari, Piran & 
Narayan 1998). For those regions where no new electrons 
are accelerated (region 5 in the stages 4 and 5, and region 
4 in the stage 5), the synchrotron spectrum is modified 
by a sharp cut-off above the cooling frequency (although 
strictly speaking the cut-off should be exponential). 

One crucial point in calculating the lightcurves is that 
the radiation from the three shocks which is emitted at 
the same fixed-frame time will reach the observer at dif- 
ferent observer times. The main reason is that the re- 
gion (4-1-5) usually moves faster (with 74) than the region 
(2-1-3) (with 73), so that for a same time interval in the 
fixed frame, dt, the emission from the sites heated by the 
shocks (B) and (C) reaches the observer (apart from the 
cosmological time dilation effect) in an interval of 

* Strictly speaking, since the region (4-|-5) trails behind the region (2 
the shell 3 ought to be taken into account. Since the shell width A3 <g; 



while the emission from the shock (A) reaches the observer 
in an interval of ^ di/",?^. In our calculations, this ef- 
fect is taken into account approximately. Defining the ob- 
server's time when the front of shell 3 reaches Rgh to be 
Tsh, we can approximate the observer's times correspond- 
ing to the regions (4+5) and (2-1-3) after Tgh to be roughly 
T4 ~ Tsh + (B-Rsh)/^lIc, andT3 - Tsh + {R-Rsh)/^llc, 
respectively**. We see that the two times are different in 
the stages 3 and 4, but are the same in stages 1, 2 and 
5. The overall effect of this correction is to "squeeze" the 
lightcurve signal from the fast moving region (4+5) during 
stages 3 and 4. 

An example of a shell collision is presented in Figure 4. 
In Fig. 4a, we plot the Lorentz factor evolution of various 
shocks and regions with respect to the blastwave radius, 
or equivalently the fixed- frame time t ~ R/c. This allows 
a clear assessment of the five stages analyzed in §3.3. Fig- 
ure 4b presents an indicative optical lightcurve, identifying 
the contributions from the three emitting sites. The final 
lightcurve is complicated, with a few features which are 
caused by different emission sites. To make more sense 
of this numerical lightcurve, it is helpful to perform some 
further analysis. The emission component from the site 
related to the shock (A) is a standard "low frequency" 
lightcurve (Sari et al. 1998), except for a gentle flatten- 
ning around 10^ s, which corresponds to the relaxation 
boosting- up of 73 in the stage 4 (Fig. 4a). The emission 
components corresponding to the shocks (B) and (C) are 
more complicated. The lightcurve slopes in various sec- 
tions depend on the (observer's) time dependence of the 
characteristic frequency, oc 'Ji'J^B' for slow-cooling 
(which is usually the case) and Uc for fast-cooling, and 
the peak spectral flux, F^y^m oc riijiB'R^. Here 7^ denotes 
the bulk Lorentz factors (73 or 74) and 7^ denotes the 
random Lorentz factors (7^, -y^, and 7^). Any change of 
the scaling law (with respect to T) of any parameter (e.g., 
7j,, 7i, rii, B', R) results in a break in the lightcurves, 
and the lightcurve slope also depends on the spectral seg- 
ment in which the observing frequency lies. In the specific 
example presented in Fig. 4b, the climbing-up segment in 
the lightcurve (C) and the flat segment in the lightcurve 
(B) correspond to the violent injection stage 3, but the 
timescale is squeezed due to the increase of 74 with respect 
to 73. The later drop-offs in both lightcurves correspond 
to the relaxation stage. Finally, all three lightcurves enter 
the post-relaxation stage at a same observer's time. There 
could be more breaks within the same stage. For example, 
the break on the lightcurve (B) in the post-relaxation stage 
is due to the transition from the (1/3) spectral segment to 
[—{p— l)/2] spectral segment. Two remarks can be made 
here. First, although in the stage 5 all the three random 
Lorentz factors decay as oc the lightcurve slopes on 

(B) and (C) in the [— (p— 1)/2] spectral segment are steeper 
than the slope on (A) in the same spectral segment. This 
arises from the different scaling law in F^^m, niore partic- 
ularly in n. For adiabatic expansion in the regions (B) 
and (C), Hi cx {V')-^ oc R-^-yi oc R-^/^ oc T-^/^ In 
the region (A), ^2 oc 73 a R~^/^ oc T~^/^. The lightcurve 
slopes for (B) and (C) then steepen by 0.75 relative to that 



-|-3), the time delay due to the emission from the region (4+5) crossing 
R, this effect may be neglected as long as (-R — Rgh) ^ ^3- 
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of (A)^. Also the whole lightcurves (B) and (C) during the 
stage 5 also tilt by 0.75, which explains why the crossing 
of Vm in the lightcurve (B) looks a little unfamiliar. Sec- 
ond, although the random Lorentz factor in the site (B) is 
much higher than those in sites (A) and (C), the flux level 
is not, mainly because the density is smaller than n2 
and during the violent collision process. 

From Fig. 4b, one can see that the injection signatures 
of a violent collision are rather complicated. There could 
be more than one signature which may be related to dif- 
ferent emission sites. This gives rise to a whole set of 
injection patterns. For example, the sharp decline of 
the site (C) emission in the relaxation stage, given an 
appropriate set of parameters, could be an explanation 
for the ROTSE prompt optical flash observation in GRB 
990123, which is conventionally interpreted as the cooling 
from the reverse-shock-heated region in an impulsive shell 
case (Meszaros & Rees 1997a, 1999; Sari & Piran 1999, 
Kobayashi 2000). The modest sticking-out of the opti- 
cal lightcurve (B) around several 10^ s (Fig. 4a) seems to 
have been observed in some GRB lightcurves, although the 
chromatic nature of this rise rules it out as a possible expla- 
nation of the achromatic bump observed in GRB 000301C, 
which can be interpreted as a gravitational micro-lensing 
event (Garnavich et al. 2000). 

The main characteristic which can be used to differenti- 
ate a post-injection model from other possible mechanisms 
leading to bumps in the light curve is the broad-band be- 
havior. Since the injection processes produce global dy- 
namics changes, they will leave imprints at all wavebands. 
In Figure 5 we present the broad-band injection lightcurves 
for (a) a violent injection; (b) a mild injection through- 
out. We can see that the mild injection case is similar to 
the Poynting-flux-dominated injection case (Fig.l), as has 
been discussed in §3.3.2. The violent injection case shows 
rather different and complex behaviors. The rising times 
in various bands are essentially the same, but the bump 
shapes are quite different, mainly because different shocks 
contribute differently in different bands. 

4. SUMMARY & DISCUSSION 

In this paper, we have analyzed the observational con- 
sequences of a GRB fireball which receives an additional 
energy injection from the central engine during the after- 
glow phase. We have discussed the cases where the injec- 
tion is either a continuous Poynting-flux-dominated wind, 
or a kinetic-energy-dominated matter shell. For the latter, 
we discuss the cases when the collision is either "mild" or 
"violent" , distinguished by whether additional shocks per- 
sist and additional emission sites arise in the outflow. We 
reach the following general conclusions: 

1. A post-burst injection energy i?inj ~ -Eimp is generally 
required to cause a noticeable signature in the afterglow 
lightcurves. This is especially so for the Poynting flux in- 
jection case and for the mild kinetic energy injection case. 
For a violent collision, a faster rear shell (i.e. 75 ^73) 
can compensate for a lower rear shell total energy (Fig. 3). 

2. For the kinetic-energy injection case, we have explic- 
itly studied the three-shell interaction process, described 
by a set of equations (eqs.[H-[||). We find that the forma- 
lin Fig. 4b, the steepening is slightly smaller than 0.75. This is 

steepened by itself. 



tion of an additional pair of shocks from the collision of 
two shells is not as common as usually thought. Depend- 
ing on the status of the rear shell when the collision occurs 
(spreading or not), a minimum relative speed between the 
two colliding shells is required, as shown in Figure 3. In 
many cases, such as expected in the internal shock model 
where slow trailing shells catch up with the fast shells when 
the latter are decelerated, the relative speed between the 
shells are usually not high enough, and the injections are 
likely to be mild. We expect that a violent injection most 
likely arises from a late injection of a high entropy shell 
from the central engine. Since the baryon loading at late 
times is expected to be smaller than that at earlier times, 
this is not unusual given a long-lived erratic central engine. 

3. If no additional emission sites form during the in- 
jection, such as the case of the Poynting-fiux-dominated 
flow, the injection signature is usally very gradual (Fig.l 
and Fig. 5b). In the case of a violent injection, the sig- 
natures are complicated and some could be very abrupt. 
In either case, the final Lorentz factor of the blastwave, as 
well as the emission flux level, are boosted up with respect 
to the case without post-injection. 

4. A prominent characteristic of the post-injection pro- 
cess is the conspiratorial variation of the fiux in all bands 
(Fig.l and Fig. 5), due to the changies in the dynamics of 
the emitting region. This provides a strong criterion for 
judging whether a bump signature is due to an injection or 
to other reasons. For example, the supernova component 
(Bloom et al. 1999; Reichart 1999; Galama et al. 2000) 
and the dust echo (Esin & Blandford 2000) interpreta- 
tions do not predict prominent changes in X-ray and radio 
lightcurves. The gravitational micro-lensing model may be 
more difficult to differentiate from a post-injection event, 
since it is essentially an achromatic feature, although it 
can have some chromatic features (Granot & Loeb 2001). 
The Poynting-flow injection case is also achromatic. The 
differentiation between the two possibilities can be only 
based on the shape of the bump. For example, the micro- 
lensing model predicts a more abrupt (shorter duration) 
feature than does the Poynting-injection. Another feature 
of micro-lensing is that the afterglow lightcurve will re- 
sume the level of no magnification when the lensing event 
ends, while in the injection case, the global energetics is 
always boosted up. Prospects of differentiation look better 
for a violent matter-injection case, in which the bumps in 
different bands have quite different shapes, although the 
rising time is essentially the same (Fig. 5a). Overall, fu- 
ture broadband afterglow data with very good temporal 
coverages plus a more detailed modelling can in principle 
differentiate the injection signature from the other possi- 
bilities. Within the injection itself, different injection sig- 
natures may provide direct information about the nature 
of the injection, and possibly, also about the nature of the 
central engine. 

The treatment presented in this paper is generic and ap- 
plicable also to more complicated injection processes. For 
example, for multi-shell-collision cases, one simply needs 
to regard the emitting gas after each injection as a new 
"impulsive" component, and a newly injected shell or wind 
can be treated in the same way. Since collisions in the in- 
ternal shock scenario also involve a similar process (but 

because the lightcurve (A) is in the (— p/2) spectral regime, which has 



12 



without the interaction with the ISM) , the treatment pre- 
sented here may be also disable in stiidying GRB gamma- 
ray Hghtcurves in the prompt phase. The shell merging is a 
very complicated process, and here we have used some ap- 
proximations leading to a semi-analytic description of the 
whole process. For example, we have assumed that the 
decelerating shell has a uniform density and speed, while 
a more accurate treatment ought to take into account the 
Blandford-McKee (1976) self-similar profiles. This approx- 
imation nonetheless provides a first-order treatment to the 
problem. Another assumption we have made is that the 
total internal energy in each region is solely due to the 
shock heating. Strictly speaking, adiabatic loss due to ex- 
pansion should be taken into account (e.g. Panaitescu et 
al. 1998). However, as long as shock heating is going 
on, the adiabatic loss is not important. Only when shock 
heating ceases (e.g., region 5 in the stages 4 and 5, and 
region 4 in the stage 5), do we consider the adiabatic loss. 
Furthermore, in all our calculations, the effects concerning 
spreading in the photon arrival time due to the curvature 
of the shells as well as their non-zero thicknesses have been 
ignored. These effects are essential in more detailed after- 
glow modeling (e.g. Granot, Piran & Sari 1999), which 
tend to smear out features on the hghtcurves. With these 
effects taken into account, we expect that the sharp fea- 
tures shown in Fig. 5a ought to be smeared out. 

Finally we compare our results with the previous rele- 
vant treatments. Recs & Mcszaros (1998), and Panaitescu 
et al. (1998) have discussed the "refreshed" injection sce- 
nario in which all the mini-shells are assumed to be ejected 
simultaneously, but with a power-law distribution of F. Al- 



though such a scenario is of interest and can be modeled 
relatively easily, it is unclear whether it is a closer rep- 
resentation of reality. The Poynting-flux-dominated in- 
jection within the pulsar central engine model has been 
discussed earlier by Dai & Lu (1998), but their indicative 
lightcurve was too abrupt. Their recent treatments (Dai 
& Lu 2001; Wang & Dai 2001) generally agree with the 
gradual evolution picture obtained in this paper. For the 
shell collision injection case, the most relevant precedent 
discussion is the one by Kumar & Piran (2000). Although 
our results generally agree with theirs on, e.g., the hydro- 
dynamic treatments, our shock forming condition (Fig. 3) 
is more demanding than theirs. In fact, for the specific ex- 
ample they presented, i.e., £"0,6 ^ ^'O.s, and 735 ^ 1.25, the 
injection is only mild according to our criterion, while they 
argue that two shocks will form and propagate into the 
two shells. This discrepancy may be due to their not hav- 
ing included the hydrodynamics of the three-shell interac- 
tion. They also presented an indicative injection lightcurve 
(their Fig.5). However, this lightcurve is more analogous 
to our mild- injection signature (Fig. 5b). As discussed in 
§3.4, a violent injection tends to generate a more compli- 
cated lightcurve with possibly more than one signature, as 
indicated in Fig. 4b and Fig. 5a. 

We are grateful to Pawan Kumar, Alin Panaitescu and 
Re'em Sari for informative correspondence and Oleg Kar- 
galtsev for discussions, to the referee for a detailed and 
helpful report, and to NASA NAG5-9192 and NAG5-9153 
for support. 
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